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2 . ABSTRACT 

Aims. Spectra, taken with the heterodyne instrument, HIFI, aboard the Herschel Space Observatory, of O-rich asymptotic giant 
branch (AGB) stars which form part of the guaranteed time key program HIFISTARS are presented. The aim of this program is to 
study the dynamical structure, mass-loss driving mechanism, and chemistry of the outflows from AGB stars as a function of chemical 
composition and initial mass. 

^ ' Methods. We used the HIFI instrument to observe nine AGB stars, mainly in the H 2 and high rotational CO lines.We investigate the 

\& , correlation between line luminosity, line ratio and mass-loss rate, line width and excitation energy. 

. Results. A total of nine different molecules, along with some of their isotopologues have been identified, covering a wide range of 

excitation temperature. Maser emission is detected in both the ortho- and para-H20 molecules. The line luminosities of ground state 
lines of ortho- and para-H 2 0, the high-J CO and NH 3 lines show a clear correlation with mass-loss rate. The line ratios of H 2 and 
NH3 relative to CO J=6-5 correlate with the mass-loss rate while ratios of higher CO lines to the 6-5 is independent of it. In most 
cases, the expansion velocity derived from the observed line width of highly excited transitions formed relatively close to the stellar 
photosphere is lower than that of lower excitation transitions, formed farther out, pointing to an accelerated outflow. In some objects, 
the vibrationally excited H 2 and SiO which probe the acceleration zone suggests the wind reaches its terminal velocity already in 
the innermost part of the envelope, i.e., the acceleration is rapid. Interestingly, for R Dor we find indications of a deceleration of the 
outflow in the region where the material has already escaped from the star. 
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1 . Introduction arising from pulsation can levitate the material and produce an 

, , r exten ded atmosphere (e.g., Bowen 1988; Hoefner & Dorfi 1992, 

Asymptotic giant branch (AGB) stars represent a late stage of [ml) Thg ^ authors showed ^ ^ models lead tQ dust 

stellar evolution when the nuclear burning in the core has ceased. formation in the atmosphere . Although the efficiency of dust for- 

The dominant factor govern ing the rest o f their evolution is the mation ifl ^ Qutflow of AGB ^ hag been tioned (see 

mass loss from the surface ( jHabing || l996fr . Presently, the initial Wolfe 2006), Hofner (2008) propo sed the wind is driven 

mechamsm(s) driving the wind which leads to dust formation in fe micron . size grains which are Fe-POor. lMattsson'eT al. ( 2008) 

O-rich stars is not fully understood. Model calculations of dy- alfjo show ^ for AGB starfj ^ a 1qw metallicit pu i sation 

namics of the photospher e of AGB stars show that Shockwaves can lead tQ an intenfje mass . loss rate Due to the larger cross sec . 

* Herschel is an ESA space observatory with science instruments pro- tions of dust grains compared to molecules, dust efficiently ab- 

vided by European-led Principal Investigator consortia and with impor- sorbs stellar radiati on and is accelerated outwa r ds, dragging the 

tant participation from NASA. gas along with it dGoldreich & Scovillg 119761 : lJusttanont et al.l 
** deceased 
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1 1994 iDecin et al.l l2006: Ra mstedt et aill2008l) . This effectively 
establishes a circumstellar envelope around the star. 

The circumstellar envelope of an AGB star is an active site 
for chemistry. The slowly expanding wind (v e ~ 10-15 km s~') 
produced by a constant, isotropic mass-loss rate, is an ideal lab- 
oratory to study physical and chemical processes during stellar 
evolution. In thermal equilibrium, the relative abundance of car- 
bon and oxygen determines the composition of the molecules 
and dust species formed. For a carbon-rich star with C/O > 1, 
carbonaceous molecules form as well as amorphous carbon and 
SiC dust after all the oxygen is locked up in the most abun- 
dant trace molecule CO. In contrast, for a star with C/O < 1, 

1. e., O-rich (M-type) AGB stars, H2O and CO are the main gas 
components, along with silicate dust. However, close to the pho- 
tosphere, shock waves due to stellar pulsations and/or the stel- 
lar radiation field can induce non-thermal e quilibrium, render - 
ing this simple picture more complicated dCherchneffl 12006). 
As an example, the canonical C-rich A GB star, IRC+10216 
has been found to exhibit H?Q emission dMelnick et al.l [200 1 ; 
lHasegawa et al.ll2006t iDecin et al.ll2010al) . Data obtained with 
Herschel- HIFI have recently shown that H2O is quite prevalent 
in C-stars dNeufeld et al.ll201 ll) . 

HIFISTARS is the guaranteed time key program aimed at 
studying circumstellar envelope s around evolved stars using the 
Heterodyne Instrument (HIFI, Ide Graauw et al.l l2010h aboard 
the Herschel Space Observatory (Pilbrattet al. 2010). The pro- 
gram aims to study stars with a large range of mass-loss rates 
(10" 7 - 10~ 4 M yr- 1 ), differing chemistry (C/0<1, C/O-l, 
C/0>1) and initial mass (low- and intermediate mass stars and 
supergiants), as well as different evolutionary (AGB to planetary 
nebula) phases. In this paper, we report all the observations done 
on O-rich AGB stars in our program. Observations and data re- 
duction are presented in Sect. 2 and the lines detected and the 
first-cut interpretation are discussed in Sect. 3. The results are 
summarised in Sect. 4. 

2. Observations 

Our HIFI observations were carried out using dual-beam switch 
mode with a throw of 3' and a slow chopping. The full bandwidth 
of HIFI of 4 GHz was utilized using the wide-beam-spectrometer 
backend. A total of 16 different frequency settings have been 
chosen to cover a number of expected strong lines of H2O and 
CO for the purpose of sampling different regions of the warm 
circumstellar envelope in objects. The data were calibrated us- 
ing the standard Herschel pipeline, HIPE, and reprocessed for 
those which had been processed with the pipeline version earlier 
than 4.0. Data were taken using two orthogonal polarizations: 
horizontal and vertical. A resulting spectrum is an average of 
these two polarizations which is then rebinned to a 1 km s _1 res- 
olution (Figs. lATl - IA. 12b . However, a number of spectra in the 
THz-band (HIFI bands 6 and 7) were affected by the ripples, 
especially in the v-polarization. These are thought to be due to 
standing waves in the hot electron bolometer (HEB) mixers. In 
these cases, we rebinned only the h-polarization spectrum. The 
spectra have been corrected for the beam efficiency, T] m b 

rjmb = rjmbfl exp(-(4no-/A) 2 ) (1) 

where <x is the surface accuracy (= 3.8 j-im), A is the 
wavelength of the transition and ?7 m b,o is a correction fac- 
tor of 0.76, except for the frequency range of 1120- 
1280 GHz (HIFI band 5) where this value is 0.66 (O lberg, 
2010, http://herschel.esac.esa.int/Docs/TechnicalNotes/ HIFI- 
Beam-Efficiencies- 1 7Nov20 1 0.pdf) 



Table 1. O-rich AGB stars in the sample. Indicated is the number 
of frequency settings observed for each object. 



name 


RA 


Dec 


obs 


D 


M 










(pc) 


(Moyr 1 ) 


IRC+lOOll 


01 06 26.0 


+ 12 35 53.0 


15 


740 


1.9E-5 (1) 


Cet 


02 19 20.8 


-02 58 39.5 


9 


107 


2.5E-7 (1) 


IK Tau 


03 53 28.9 


+ 11 24 21.7 


15 


260 


4.5E-6 (1) 


R Dor 


04 36 45.6 


-62 04 37.8 


9 


59 


2.0E-7 (2) 


TXCam 


05 00 50.4 


+56 10 52.6 


9 


380 


6.5E-6 (1) 


WHya 


13 49 02.0 


-28 22 03.5 


16 


77 


7.8E-8 (1) 


AFGL 5379 


17 44 24.0 


-31 55 35.5 


8 


580 


2.0E-4 (3) 


OH 26.5+0.6 


18 37 32.5 


-05 23 59.2 


8 


1370 


2.6E-4 (3) 


RCas 


23 58 24.9 


+51 23 19.7 


9 


106 


4.0E-7 (1) 



Referenc es. The quoted dist ance and mass-loss rates are taken 
from (1) | de Beck et all d201Ch : (2) iMaercker et all i2009h : and (3) 
Uusttanont et alj ( 120061) 



In three objects, IRC+ 10011, IK Tau, and W Hya, we ob- 
serve a large number of frequency settings (15, 15 and 16, re- 
spectively), with the intention of using them as templates for 
other objects with high, intermediate and low mass-loss rates, 
respectively. In order to study a larger sample of stars, we also 
observed 6 further objects with fewer frequency settings selected 
from the original settings (TableQ} which include many diagnos- 
tic lines (see Table lA.ll . 

In order to estimate the integrated line intensity, a Gaussian 
fit was performed in IDL. Some lines are flat-topped and hence 
the line intensity was estimated using a rectangular (rhombic) 
profile. Due to the observing mode used, it is not possible to 
recover the true continuum level of the spectra. However, the es- 
timated line intensities are not affected by this since we subtract 
the continuum prior to calculating them. Our error estimates re- 
flect only the noise in the baseline. Other uncertainties such as 
the assumption that a line is Gaussian are not included. However, 
independent measurements of the line intensity by summing the 
area under the line agree within twice the estimated uncertain- 
ties listed in Table lA.ll The absolute flux calibration error is ex- 
pected to be ~ 30% in the HIFI bands 6 and 7 and less (~ 15%) 
at lower frequencies. 



3. Molecular lines 

We identify a total of 9 different molecular species as well as 
associated isotopologues (Table IA.U . Most of the lines are in 
the ground vibrational state. However, for H2O and SiO, vibra- 
tionally excited lines are also seen indicating that the lines orig- 
inate from the hotter part of the envelope (Sect. 13.31 ). Many of 
the settings observe d in the star IK Tau have been presented by 
IDecin et alJ ((2010b) and are included here for completeness of 
our sample, along with new observations. 

It has been postulated that H2O is one of the main molecu- 
lar coolants in the circumstellar envelopes of O-rich AGB stars, 
along with CO (Goldreich & Sco villell976l) . This was confirmed 
by observations by the Infrared Space Observatory (ISO) of 
these stars wh ich show numerous strong lines of b oth these 
molecules (e.g. lBarlow et alJll996LlNeufeld et al.ll 19961) . We de- 
tected a total of 23 H2O lines in all the three main isotopologues 
covering the excitation temperature from 30 to 2350 K. These 
lines probe the full extent of the circumstellar envelope and due 
to HIFI resolution, the lines are well resolved, enabling us to 
study the dynamics of the dust-driven wind. It should be noted 
that the H2O 1 i,o _ lo,i vibrationally excited line at 658.0GHz is 
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Fig. 1. Integrated line intensities for all stars in the sample. CO 
: black open circles; 13 CO red filled circles; H^O : open squares; 
H^O : red filled squares; Hl s O : green filled squares; SiO : stars; 
SO : open triangles; SO2 : x; HCN : filled triangles. 



Fig. 2. A plot of the line luminosity (squares : ortho-H^O 1 10- 
lo,i, circles : para-H^O 1 i,i-Oo,o, stars : CO J=6-5, x : CO J=10- 
9, triangles : NH3 Iq-Oo) as a function of the published mass-loss 
rate show a strong con-elation between the two parameters. 



likely to be a maser, as well as the 620.7 GHz (53,2-44,0 line (see 
Sect. 13. 2\ . This latte r line has previous ly been reported in the su- 
pergiant VY CMa (lHarwit et alJ l2010). As for CO, we detected 
a total of 8 transitions in three isotopologues. 

Other molecules are also identified, such as NH3, SiO, HCN, 
SO, SO2 and OH. These molecules have already been detected 



in HI FI spectra of AGB a nd post- AG B stars (Buiarrabal et al 
2010s lDecinetalJl2010tJ ; iMenten et all 120101; I Justtanont etal 
20ld ISchoier et al.ll201 lb . In AFGL 5379, a line is detected at 



1 196.010GHz which can be attributed to the H 2 S (3,1,2)-(2,2,1) 
transition. This line is not seen in the other stars in our sam- 
ple, however it is present in the supergiant VY CMa (Alcolea 
2011, in preparation). A line is also seen in the spectrum of 
OH 26.5+0.6 at 1 1 14.431 GHz which corresponds to the vibra- 
tionally excited SiS J=62-61 transition. Although, with the high 
upper energy level of 1918.3836cm -1 , we classify this line as 
unlikely. 

Since HIFI employs a double-side band mode, there is 
a possibility of ambiguity in classifying a line. The line at 
1112.833 GHz in the lower side-band due to 29 SiO J=26-25 co- 
incides with the H2O 116,6-123,9 at 1 101.130GHz. Considering 
the upper energy levels and the expected line strengths of the 
two, we conclude that the line is likely due to 29 SiO. 

3.1. Line intensities 

All objects show strong emission due to H2O and CO molecules 
(Fig. Q). The OH/IR stars, OH 26.5+0.6, AFGL 5379 and 
IRC+10011, show very strong emission of both ground state 
ortho- and para-H20 lines compared to that of CO J=6-5. 
Despite the high mass-loss rates in these OH/IR stars, the CO 
J= 16- 15 is very weak or not detected. This may be partly due 
to the attenuation of stellar light by dust, preventing the excita- 
tion of this line. In o Cet, the CO lines are much stronger than 



other molecular species. This is lik ely due to the white- dwarf 
companion with its hard X-ray flux (lO'Dwveret al.ll2003l) pho- 
todissociating H2O molecules in the circumstellar envelope of 
the primary star. In all the stars, we detect both 12 CO and 13 CO. 
The line ratios of the corresponding transitions vary from 1.5 
in stars with high mass-loss rate to up to 10 in stars with lower 
mass-loss rates. 

We detected both the ortho- and para-H20 as well as the iso- 
topologues and vibrationally excited lines. A discussion on iso- 
topic ratios of detected lines is presented by Silva et al. (201 1, in 
preparation). A detailed modelling of radiative transfer of stars 
with low mass-loss rates is in progress (Maercker et al. 201 1, in 
preparation). The highest excitation lines seen in our spectra are 
the vibrationally excited 2 12 - lot and the vibrationally excited 
ground state lines of both ortho- and para-H20 with excitation 
energies in excess of 2000 K (Table |A~TT >. These lines are likely 
radiatively pumped. 

In Fig. |2 we plot the line luminosity as a function of mass- 
loss rate for CO J=6-5, 10-9, and the ground-state lines of 
both ortho- and para-H20, as well as for NH3 1q-0q. The dis - 
tances and mas s -loss r at es are taken from Ide Beck et al.1 ([2010); 
iMaercker etal] (120091) ; Uusttanont etal] J2OO6I) (see Table [Q. 
There is a strong positive correlation between individual H2O 
and CO line luminosities, i.e., the slopes for the mass-loss rate 
up to ~ 10~ 5 M Q yr~' are betwen 0.7-0.8. There is one excep- 
tion regarding the H2O luminosity in o Cet which is much lower 
than for the CO. This is likely due to H2O being photodisso- 
ciated by the binary companion. The flattening off of this rela- 
tion at M ~ 3 x 10 ~ 5 M yr~' is consistent with relations by 
Ide Becketal](l2010l) . 

Another molecule readily detecte d in circumstellar en- 
velopes of AGB stars is SiO (e.g., Buiarr abal et al.l 1.1994; 
Olof sson et al. I IT998I Here, we identified three silicon isotopes. 
The vibrationally excited lines of high- J transitions are also seen. 
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Fig. 3. A plot of the line ratios relative to the integrated inten- 
sity of the CO J=6-5 as a function of mass-loss rates. In the 
lower panel : for ortho-FbO 1 1,0- lo,i (squares), para-H20 1 1,1 - 
Oo,o (circles), NH3 lo-Oo (triangles) scaled up by a factor of 5. 
The the middle panel shows CO J=10-9 (x), CO J= 16- 15 (stars) 
and in the upper panel : SO ( 1 3 13- 12.12) 559.3 and SO (13 14- 
12 13 ) 560.2 GHz (triangles) and S0 2 (21 6 ,i 6 -21 5 ,i 7 ) 558.4 and 
S0 2 (37i,37-36 ,36) 659.4GHz (squares). 



SO has been observed in a number of AGB stars covering a 
wide range in mass-loss rat es. The estimated abundance is a few 
1(T 7 up to 1(T 6 M yr _1 dOmont etalj|1993l iBuiarrabal etall 
[1991 . However, we only observed SO and SO2 lines in our 
HIFI frequency range in stars with low mass-loss rates. From 
ISO observations, SO2 ab sorption is detected tow ards stars with 
optically thin envelopes dYamamura et al.l fl999), while this is 
not seen in optically thick shells. SO is thought to be formed 
mainl y via S + CO and to a lesser extent via S + OH dCherchneffl 
2006)' and is expected to be present throughout the envelope. It 
follows then that the formation of SO2 is via SO + OH. The 
best example for SO and SO2 can be seen in our spectrum of R 
Dor (Fig. IA.6l l. In this object, the line fluxes of the SO lines are 
about an order of magnitude larger than that of SO2. In the other 
stars where both molecules are observed, the flux ratios are less 
than 5. The only other S-bearing molecule seen in our spectra is 
H 2 S in AFGL 5379 at 1 196.0 GHz. This line is absent in all the 
other stars in our sample, compare to a more ubiquitous pres- 
ence of this molecu le observed from the ground in OH/IR stars 
dOmont etal.lll993l) . 

We detected the OH triplet line at 1834.7 GHz, which is very 
strong in the two extreme OH/IR stars, in agreement with the fact 
that these stars emit strongly in the OH 1612 MHz maser (e.g., 
ISevenster et al.ll200l1) . The spectrum of IRC+1001 1 is contami- 
nated by ripples (standing waves in the HEB mixers) preventing 
us from confirming the presence of the OH line. 

To study the effects of mass-loss rate on the excitation, we 
plot line intensity ratios of CO J=16-15 and J=10-9 relative to 
the J=6-5 (Fig. [3] middle panel). These ratios appear to be con- 
stant over four orders of magnitude of mass-loss rates, indicating 



Fig. 4. The lower panels show a comparison of the normalized 
line profile of the 620.7 GHz H2O 53,2-444 line (solid line) with 
the ground state ortho-H20 li,o-lo,i li ne at 556.9 GHz (dotted 
red line) and the CO J=6-5 (dashed blue line). The upper panels 
show the 970.3 GHz H 2 5 2 ,4-4 3 ,i (solid) with same the H 2 
ground state and CO lines. 



that the density is higher than critical density of these transi- 
tions, i.e., the lines are thermalized in the emitting region. The 
line ratios of the ground state lines of both ortho- and para-H20 
with CO, on the other hand, show an increasing trend as a func- 
tion of mass-loss rate with the para-H20 having a consistently 
higher ratio than that of the ortho-H20. This same increasing 
trend can be seen between the NH3 lo-Oo and mass-loss rates, 
indicative of similar excitation conditions for both molecules. 
The top panel of Fig. |3] shows the SO/CO and SO2/CO ratios. 
Here, a decreasing trend for the line ratio with the mass-loss 
rate can be seen. From the HIFI observations, both SO and SO2 
are not detected in stars with a mass-loss rate higher than 10~ 5 
M yr~' . It i s therefore not possibl e to confirm the tight correla- 
tion seen bv lOlofsson et all d 19981) b e tween CO (J=l-0) and SO 
(Jk=32-2i) line fluxes. lOmont et al.1 d 19931) . however, detected 
low excitation lines of SO and S0 2 (T ex < 55 K) in 14 OH/IR 
stars, including three stars in our sample (IRC+1001 1, IK Tau 
and OH 26.5+0.6). In optically thick envelopes, the high exci- 
tation levels of the sulphur-bearing molecules in the HIFI range 
(T ex > 200 K) are simply not sufficiently excited. 

3.2. H 2 masers 

Maser emission has been observed toward a large number of 
evolved stars. Its signature includes an anomalously strong line 
intensity and a narrow line width due to maser amplification. In 
some instances, a maser can be seen in a narrow peak in the blue 
wing. 

In three objects (W Hya, IK Tau, and IRC+1001 1), we ob- 
serve the 620.7 GHz tra nsition of ortho-H^O 5^.7-4 4. 1 which is 
predicted to be masing (Neu feld & Melnicklll991l) . The maser 
line is brighter compared to the thermal excited lines and tends 
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Fig. 5. A comparison of the normalized line profile of the 
658.0 GHz H 2 line (solid line) with the ground state ortho-H 2 
line at 556.9 GHz (dotted red line) and the CO J=6-5 (dashed 
blue line). 



to have a narrower line profile as the region of coherent veloc- 
ity required to produce maser emission is generally small. In W 
Hya, the maser line is narrow compared to the thermally excited 
lines of H 2 1 10 - lo,i and CO J=6-5, indicating it comes from 
a region close to the star. In IK Tau, the line exhibits a double- 
peak profile which has a line width slightly narrower than that 
of the 556.9 GHz ground-state ortho-H 2 and the CO line. In 
IRC+1001 1, however, the line is of the same width as the ground 
state H 2 and CO J=6-5 lines, and is thought to arise from ther- 
mal emission, but it also has a strong blue component due to the 
masing action (see Fig. |4]lower panels). This is in line with the 
expectation t hat the 620.7 GHz line is masing at high (~ 10 3 K) 
temperature (Neufeld & Mel nickl 19911) . The maser emission ap- 
pears to arise from the blue-shi fted part of the spectrum, i .e., in 
the front side of the envelope. Neufeld & Mel nickl (Il99ll) also 
predicted a maser line at 970.3 GHz from para-H 2 (5 2j 4-43 1) 
which we observed in W Hya and IK Tau with similar line widths 
and profiles as for the 620.7 GHz maser. In IRC+1001 1, we did 
not detect a maser peak (Fig. |4] upper panels), and the line is 
about the same brightness as for the ground transition line at 
556.9 GHz (Table lATTT i. hence the maser action appears to be 
quenched in this object. 

In one of the frequency settings, we detected the vibra- 
tionally excited H 2 li,o-lo,i li ne at 658.0GHz which is likely 
masing (Fig. 01. In all the objects with low and intermediate 
mass-loss rates, the maser line is markedly narrow compared to 
the 556.9 GHz line except for o Cet, where both lines have com- 
parable widths. This indicates, again, that the maser line in the 
majority of the stars originates close to the star where the wind 
has not yet reached the terminal velocity. The intrinsic line inten- 
sity supports the fact that the line is masing, i.e., the line intensity 
of the vibrationally excited line is very bright, compared with 
the ground state line at 556.9 GHz. However, in extreme OH/IR 
stars (OH 26.5+0.6 and AFGL 5379), the widths of the ground 



and vibrationally excited lines are comparable and the flux of 
the vibrationally excited line is small relative to the ground state 
line. It is unlikely then in these two objects that the line is mas- 
ing. This is probably due to the high density quenching maser 
action in the inner part of the envelope. 

We also observe the vibrationally excited line of para-H 2 
li,i-Oo,o at 1205.8 GHz. Unlike the vibrationally excited 1 i,o _ lo,i 
ortho-H 2 0, this line does not appear to be masing when compar- 
ing with the same transition in the ground vibrational state at 
1 1 13.3 GHz as the vibrationally excited line is much weaker. 

3.3. Expansion velocity of circumstellar envelops 

Once dust grains condense in the outflow of AGB stars, the 
wind acceleration mechanism is thought to be due to radiation 
pressure on the dust grains, which drags the gas as they move 
away from the cen tral star (e.g., iGoldreich & Scovilld 1 1976b 
Habingetal. 1994). Lines with a high excitation temperature 
are expected to be emitted in regions close to the star, where 
the wind has not yet reached its terminal velocity. Consequently, 
these lines would be observed to be narrower than lines with 
lower exc itation temperatures, e mitted from farther out in the 
envelope. lJusttanont et al.l ([1994) demonstrated that in a radia- 
tive transfer calculation for successively higher-J CO lines, the 
warmer regions of the emitting zone are being probed which 
have smaller line widths than for the lowest CO transition. 
However, observations of maser lines which probe regions close 
to the central star revealed that the wind acceleration is not a s 
fast as predicted by the dust-drag wind (Chap man et al.| [l994). 
In Fig. [6] we plot the expansion velocities for the lines detected 
in our spectra, which probe the warm region as well as the region 
close to the central star. Here, the observed expansion velocity is 
a measure of half the width at the baseline level. The uncertainty 
of the expansion velocity is ~ 15-30%. In most cases, there is 
a trend that lower excitation lines have wider observed line pro- 
files, indicating they come from regions where the wind is close 
to the previously determined expansion velocity from ground- 
based observations of low-J CO lines. This indicates that there 
is a velocity gradient in the circumstellar outflow. One caveat 
is that for H 2 lines, the profile can be heavily absorbed in the 
blue wing, leading to slightly lower measured expansion veloc- 
ities. This can be seen in Fig. [6] which shows that the measured 
expansion velocity of CO J=6-5 is larger than that for the ground 
state of both ortho- and para-H 2 0. We observe isotopic lines of 
H 2 0, and a trend can be seen that v(H 2 e O) > v(H 2 7 0, H 2 s O). 
This may be due to the fact that the less abundant species are 
more prone to photodissociation from the interstellar radiation 
hence we see molecules close to the central star where the gas 
has not yet reached the terminal velocity. 

For all stars, 3 CO lines have been targeted : J=6-5, 10-9 
and 16-15. The widths of these lines show the trend of increas- 
ing velocity for the lower excitation line, in agreement with the 
acceleration. However, the J=16-15 line is not seen towards the 
two extreme OH/IR stars, OH 26.5+0.6 and AFGL 5379, possi- 
bly due to the dust attenuation of stellar radiation, preventing the 
excitation of this line. We readily detected the 13 CO lines in our 
frequency settings as well and in most cases, the line widths of 
these two isotopes are comparable. 

The line widths of SO and S0 2 are narrower than those from 
CO and H 2 0. Due to the low abundance, the lines are not thought 
to be optically thick, hence the reason why the lines are nar- 
rower is likely because they originate from the inner region of 
the envelope. It can be that these molecules are also destroyed 
closer in the central star by interstellar UV radiation, compared 
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Fig. 6. Observed expansion velocity for molec- 
ular lines in our sample stars. 



to the more robust CO molecules. In some stars with a relatively 
low mass-loss rate, the measured velocities of SiO J=16-15 and 
2y SiO J=13-12 lines also are smaller than those from CO. 

A number of stars show a clear trend of the high exci- 
tation lines having a smaller line width than low excitation 
ones (Fig. [§}, consistent with a region where the molecules are 
then accelerated towards the terminal velocity observed in the 
ground-based low-J CO observations. In order to qualitatively 
show that the wind is acce lerated, we overplot a functi on similar 
to the dust-driven wind o f lLamers & Cassinelli ( 1999) but with 
a dependence on the energy rather than the radius : 



Table 2. Parameters describing the velocity profile for stars 
which show wind acceleration (Fig.[6]l. 



v = vo + (v e -vo)[l.-Eo/Ef 



(2) 



where vrj, v e and Eq are arbitrary initial and expansion veloci- 
ties in km s _1 and energy in K, respectively, and /3 is a veloc- 
ity exponent (Table 0. The parameters are chosen to follow the 
12 CO transitions as a the main tracer for the expansion veloc- 
ity. In a few cases, however, no clear jump in velocity between 
high and low excitation lines is seen, and a straight line can bet- 
ter describe the velocity field (right panels in Fig. [6). The fitted 
lines are least-square-fits to all the data points. In general, de- 
spite the smaller observed expansion velocity of SO, SO2, SiO 
and less abundant isotopologues of H2O and CO, the measured 
velocities of these species follow the overall trend of increasing 
velocity with decreasing excitation energy, i.e., increasing radius 
from the central star. 

In the case of R Dor, however, the highly excited lines of 
both H2O and SiO indicate a larger expansion velocity in the in- 
ner part of the envelope than for the lower J-transitions of CO. 
This is counter-intuitive compared to what is expected of an ac- 
celerating dust-driven wind and may pose a challenge to detailed 
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radiative transfer modelling of the stellar wind of R Dor. Also in 
AFGL 5379 and OH 26.5+0.6, the expansion velocity of var- 
ious lines suggests that the wind reaches the terminal velocity 
already in the innermost part, i.e., the acceleration zone is very 
small since the line widths are similar regardless of the excita- 
tion. For W Hya, the line widths increase linearly with the de- 
creasing excitation energy. These stars demonstrate that the gas- 
dust interaction does not seem to follow a simple momentum 
transfer previously used to describe the wind. To fully under- 
stand the observed line widths, a full radiative transfer must be 
performed including a velocity field and temperature profile. 



4. Summary 

We present the full dataset on O-rich AGB stars as part of 
the guaranteed time key program HIFISTARS to study the late 
stages of stellar evolution both kinematically and dynamically. 
We detect emission of 9 molecular species and their most com- 
mon isotopologues. We find a trend of increasing line luminosity 
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of H2O, CO and NH3 with the mass-loss rate (Fig. [2]) which ta- 
pers off for stars with a high mass-loss rate. Interestingly, the line 
ratios of high-J CO relative to CO J=6-5 are independent of the 
mass-loss rate (Fig. [3}. This implies that the excitation of high-J 
CO is independent of the density, i.e., the lines are fhermalized. 
Other species such as H2O and NH3 show a positive correlation 
with the mass-loss rate. 

From the line brightness and shape, we conclude that the 
H2O 620.7 GHz is masing in all three objects observed. The 
970.3 GHz and the 658.0 GHz lines are masing in objects with 
relatively low (< 10~ 5 M Q yr -1 ) mass-loss rates. These maser 
lines are generally significantly narrower than the thermal H2O 
lines. 

Since the lines are well resolved, we can use the observed 
line width to derive the expansion velocity for each transition. 
In general, highly excited lines which originate close to the star 
have smaller expansion velocities compared to the lower exci- 
tation lines (Fig. [6}. This is in agreement with the gas being 
dragged by dust grains as they are accelerated outward due to 
stellar radiation. However, from our observations, the accelera- 
tion zone differs in our sample, with stars with a high mass-loss 
rate showing high-excitation lines having velocity close to the 
terminal velocity while stars with a low mass-loss rate having a 
distinct low velocity region (E u < 10 3 K). In R Dor, the highly 
excited lines show smaller observed line widths than the low ex- 
cited lines. For this object, the gas appears to be decelerating as 
it moves away from the central star. A caveat is the optical depth 
which can affect the line profile. In order to calculate this effect, 
a full radiative transfer has to be performed in individual objects. 
This work has started and will be presented in the forth coming 
papers. 
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Appendix A: Spectra 

In this section, we present the observed spectra of O-rich AGB 
stars observed in the HIFISTARS guaranteed time program, 
along with the measured line intensities. 
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Fig. A.10. HIFI spectra of AFGL 5379. 
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Fig. A.11. HIFI spectra of OH 26.5+0.6. 
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Table A.l. Calculated line intensities (K km s ') of the observed lines. The setting numbers are internal HIFISTARS LO frequency settings (see also Figs IA.1HA. 121 . The 
marks the settings which have not been observed and 'x' indicates that the line is not detected above the noise level. 
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